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Abstract 

Model simulations of solar irradiance reaching the 
Earth's surface during a solar eclipse constitute a useful 
tool for studying the impact of this phenomenon on the 
radiance propagation through the atmosphere. A sim- 
ple approach to extend the use of an algorithm already 
adopted for evaluating the variations in the extraterres- 
trial solar radiation during a total eclipse is proposed 
for a partial eclipse case. The application is based on 
the assessment of the distance between the apparent so- 
lar and lunar disk centers on the celestial hemisphere, 
using the local circumstances and the ratio between the 
Sun and Moon radii as input parameters. It was found 
that during the eclipse of March 29, 2006, the present 
approach led to an estimate of the surface UV solar 
irradiance trend differing by no more than ±5% from 
the corresponding trend observed at Bologna (Italy). 

1 Introduction 

A solar eclipse is a good opportunity for study- 
ing the interaction between solar radiation and the 
Earth's atmosphere, [TJ HJ El HI E], and for test- 
ing radiative transfer models [6]. Together with 
ground-based measurements, modeling of the in- 
coming solar irradiance at the surface during an 
eclipse episode can provide useful information for 
those investigating such interaction processes [5] . 

Mollman and Vollmer [7] proposed an approach 
aimed to evaluate the broadband solar irradiance 
during a solar eclipse, in which the limb darken- 



ing effect was neglected. An accurate procedure 
for calculating the spectral variations in the extra- 
terrestrial solar radiation occurring during a total 
eclipse was developed by Koepke et al. [5] ■ The ap- 
proach takes into account the limb darkening and 
determines the radiance changes in terms of nor- 
malized spectral irradiance I norm ,\ expressed as 
the ratio between the radiation coming from the 
part of the solar disk uncovered by the Moon and 
the radiation emitted by the entire solar disk. The 
parameter I nor m,\ characterizes the variations in 
solar irradiance at the top of the atmosphere during 
an eclipse episode and, therefore, it can be used as a 
correction factor for evaluating the extra-terrestrial 
irradiance in the radiative transfer models. The 
main parameter describing the eclipse geometry in 
the method is the distance between the apparent 
solar and lunar disks, which was evaluated in the 
total eclipse case. The present study aims to ex- 
tend the approach developed by Koepke et al. [5] 
to a partial eclipse case, assessing the length of 
segment defined by the Sun and Moon centers as 
a function of the ratio between the corresponding 
radii and local eclipse circumstances. 

2 Assessment of the solar-moon ap- 
parent distance during a partial 
solar eclipse 

In case of total eclipse, the apparent path of the lu- 
nar disk center on the celestial hemisphere, crosses 
the corresponding center of the solar disk, and pa- 
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Fi gUre 1: Apparent movement of the Moon disk of radius 
Rm, labeled with apex (') for the first phase of the eclipse and 
apex (") for the second, with respect to the Sun disk of radius 
Rs , during a partial solar eclipse episode. The dashed circles 
represent the position of the Moon at times relative to (from 
right to left): (i) first contact, (ii) maximum eclipse, and (iii) 
fourth contact. The distance X between solar and lunar disk 
centers is similarly labeled with apices ( ) and (") for two phases 
of the eclipse, respectively. D is the distance between the Sun 
and Moon disk centers at the maximum eclipse time £ m ; M is 
the part of the solar disk radius covered by the Moon at time 
t m ; and S (and S at time t) are the apparent paths of the 
Moon disk center during the eclipse, labeled with apices (') and 
("), respectively, as for the above-mentioned parameters. 



rameter I„orm,\ can be determined [5] as a func- 
tion of the distance X between the Sun and Moon 
centers, which is assumed to follow a linear time- 
dependent trend for a constant velocity v of the 
Moon with respect to the Sun. However, in the 
case of a partial solar eclipse, the Moon disk cen- 
tre does not cross the corresponding solar center 
[S], as is shown in Fig. [T]and, therefore, parame- 
ter X does not vary linearly in time. The graph 
represents the relative movements of the apparent 
solar and lunar disks of radii Rs and Rm, respec- 
tively, during an eclipse episode starting at time 
t Q , reaching its maximum at time t m and finish- 
ing at time t e . All the distances shown in Fig. Q] 
are considered to be normalized to parameter Rs, 
so that the apparent Sun disk radius is considered 
as unit (Rs = 1). In addition, the solar disk has 
been assumed to have a fixed position. Due to the 
specific geometry of the problem, the time-interval 
between t and t m is usually slightly different from 
the corresponding time- interval between t m and t e . 
This is why the eclipse episode is arbitrarily subdi- 
vided into two parts, the former of which is taken 
between times t Q and t m and labeled with apex ('), 
and the latter between times t m and t e with apex 
("). According to Fig. [TJ parameter X' can there- 
fore be expressed as: 



X' = ^D 2 + S' 2 . 



(1) 



Parameter D is the distance between the centers 
of apparent Sun and Moon disks at the maximum 
eclipse time and can be determined in terms of the 
following equation, 



D =R S + R° 



M 



M 



1 + R° M -M, (2) 



where R M is the Moon radius at time t m . Seg- 
ment M, which represents the part of solar radius 
Rs covered by the Moon at the maximum eclipse 
time t m , can be expressed in terms of the eclipse 
magnitude M g for the considered site in the follow- 
ing form: 



M = 2R s M q 



2M„ 



(3) 



At the eclipse start time t , parameter X' Q is equal 
to Rs + R'm an d, hence, during the period [t ,t m ], 
the Moon disk centre describes the path S' shown 
in Fig. [TJ whose length is given by 



S' = J{Rs + R' M Y-D 2 



= ^(l + R' M r-D*. (4) 

Assuming nearly constant velocity v' of the Moon 
with respect to the Sun during the first phase of 
the eclipse, it can be stated that 



(5) 



Therefore, parameter S' can be defined at each 
time t € [toj^m] m the explicit form, 



S' = S' - v' (t - i ) 



(6) 



Thus, the distance X' can be defined using Eqs ([T]) 
- ([6]). Similarly, the corresponding parameter S" 
can be expressed as equal to the difference: 



S" = Sq - v'o(t e - t) 



(7) 



where t € {t m ,t e ] and parameters Sq and v'q are 
defined according to the following pair of equations: 



and 



(1 + R' M ) 2 -D 2 



on 



(8) 
(9) 
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Table 1: Parameters of the solar eclipse of 29 March 2006 
estimated by Espcnak and Anderson 1121 for Bologna, Italy. 
The values of ratio Rm/Rs arc taken from physical ephemcris 
of the umbra! shadow (Table 4 in |12| ). 



Contact 


Time 
(UTC) 


Rm/Rs 
{Rs = 1) 


Magnitude 

Mg 


I 


09:33:01 (t D ) 


1.0500 (R' M ) 




Maximum 


10:38:03 (f m ) 


1.0509 (R" M ) 


0.532 


eclipse 








IV 


11:44:11 (t e ) 


1.0375 {R'lf) 





assuming again that velocity v'q is nearly constant. 
Finally, the distance X" between the apparent so- 
lar and lunar disk centers during the second phase 
of the eclipse is defined according to Eq ([T]), as 

X" = VD 2 + S" 2 . (10) 

Equations ([1]) - (fTU|) provide an algorithm suit- 
able for determining the time-patterns of parame- 
ter X, denoted as X' during the first half of the 
event and as X" during the second half, know- 
ing the local circumstance times, eclipse magnitude 
M g and behavior of the ratio Rm/Rs during the 
eclipse. In fact, the present procedure takes into ac- 
count the ratio Rm /Rs only at the start time t Q of 
the event, at the maximum eclipse time t m , and at 
the final time t e , which correspond to parameters 
R' M , R° m and R M , respectively, assuming Rs = 1 
as established above. Since the variations in ra- 
tio Rm/Rs usually do not exceed 1-2%, the error 
made in estimating parameter X when considering 
such a ratio at three fixed times only, will be around 
1%. In case of total eclipse, parameter D is equal 
to zero and the above approach gives a linear trend 
of distance X, as assumed by Koepke et al. [5] It 
is also worth noting that Koepke ct al. [8] normal- 
ized the parameter X by Rm + Rs in their proce- 
dure, while such a step was avoided in the present 
procedure for sake of simplicity. The above algo- 
rithm is similar to that proposed by Mollman and 
Vollmer [7] to evaluate the illuminance during a so- 
lar eclipse. The main difference is the subdivision 
of the event into two periods, before the maximum 
phase and after that, while Mollman and Vollmer 
[7J assumed a symmetry of the eclipse episode evo- 
lutionary features with respect to the maximum. 
On the other hand they assessed the broad-band 
irradiance and, hence, the limb darkening was not 
taken into account. Conversely, the Koepke et al. 
[5] method considers this effect: based on these 
concepts, the present algorithm is suitable to esti- 
mate the spectral variations of surface solar irradi- 
ance during a partial solar eclipse. 




9:30 10:00 10:30 11:00 11:30 12:30 



TIME (UTC) 

Figure 2: Time-patterns of parameters X and Inorm. 325 dur- 
ing the partial solar eclipse of 29 March 2006 evaluated for 
Bologna (Italy): X is the distance between the Moon and Sun 
disk centers; I norrTlj 325 is the normalized solar irradiance at 
the 325 nm wavelength, given by the ratio between the radia- 
tion coming from the uncovered part of the solar disk and that 
coming from the whole disk. The contact times and maximum 
eclipse time (see Table 1) arc indicated by vertical dotted lines. 

3 Validation 

To illustrate the outputs given by the present 
algorithm, the global solar UV irradiance time- 
patterns observed during the partial solar eclipse 
of March 29, 2006 at Bologna, Italy (44°31'N, 
11°20'E), were evaluated and compared with the 
ground-based measurements [10] performed with 
the narrow- band filter radiometer UV-RAD [IT] . 
TablcJT] presents the input parameters, as provided 
by Espenac and Anderson fT3|. Figure [2] shows 
the time-patterns of normalized solar irradiance 
,325 at the 325 nm wavelength, as calculated 
following the Koepke et al. [8] procedure, together 
with those of distance X evaluated through Eqs 
(HJ) - ([TU]) . Parameter I n0 rm,325 defines the relative 
variations of the radiance entering the atmosphere 
at the 325 nm wavelength, occurring at Bologna 
during the event. It is suitable for use as a correc- 
tion factor for extra-terrestrial solar radiation in a 
radiative transfer model used to evaluate the sur- 
face solar irradiance. To perform such an evalua- 
tion, the widely used TUV model [12] was applied 
for the 325 nm spectral component of solar irra- 
diance, keeping constant the daily mean value of 
columnar ozone amount equal to 390 DU, as de- 
termined by UV-RAD on March 29, 2006. Since 
the eclipse day was characterized by stable clear- 
sky conditions, cloud effects were not considered in 
the present assessment. On the basis of ground- 
level visibility evaluations made on the eclipse day 
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|10j . an aerosol volume extinction coefficient was 
assumed at the 550 nm wavelength varying between 
less than 0.04 and 0.05 km -1 . Consequently, as- 
suming a particle scale height close to 1.25 km ac- 
cording to the Penndorf [14| evaluations (confirmed 
by the estimates of 1.0 - 1.5 km made by Tomasi 
[15] in the Po Valley on early spring days), aerosol 
optical depth (AOD) at visible wavelengths can be 
evaluated to range between 0.04 and 0.07, for me- 
teorological and atmospheric transparency condi- 
tions similar to those of the eclipse day. Thus, 
an AOD value equal to 0.05 was given as input 
to the code, together with the Elterman [16] verti- 
cal distribution of the volume extinction coefficient 
of aerosol particles. The winter Mid-latitude tem- 
perature profile determined by Anderson et al. [T7] 
was used in the code calculations for surface albedo 
equal to 0.08, leading to the best agreement be- 
tween evaluated and measured irradiances during 
the part of the day in which perturbations due to 
the eclipse were not observed. 
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Figure 3 1 Upper part: time-patterns of surface UV solar ir- 
radiancc at the 325 nm wavelength, measured during the par- 
tial solar eclipse of 26 March 2006 at Bologna (open circles). 
The thick grey curve represents the corresponding model eval- 
uations. Lower part: time-patterns of the ratio between UV- 
RAD measured and model-evaluated irradiances, respectively. 
The contact times and the maximum eclipse time (sec Table fit 
are indicated by vertical dotted lines. 



In the upper part of Fig. [3] the results of the 



simulation are presented, showing a comparison be- 
tween the theoretical and observed variations in so- 
lar irradiance. The lower part of Fig. [3]prcsents the 
time-patterns of the ratio between measured and 
calculated irradiance values, providing evidence of 
the good agreement between the two variables. Ex- 
cept for the times of maximum eclipse and last con- 
tact, when discrepancies of about -5% were found, 
all the other cases exhibit similar differences within 
±3%. The comparatively higher deviations be- 
tween evaluated and measured values at times t m 
and t e can be reasonably attributed to the presence 
of sparse cumulus congestus clouds moving around 
the Sun on 29 March 2006, causing, together with 
the eclipse event, appreciable changes in the sur- 
face solar irradiance that cannot be represented 
with good precision by using a monodimensional 
radiative transfer model. 

4 Conclusions 

A simple algorithm has been proposed suitable for 
extending the method of Koepke et al. [5], devel- 
oped to correct the trend of extra-terrestrial solar 
radiation during a total solar eclipse to a partial 
eclipse case. To do this, the evolution of the dis- 
tance between the apparent Moon and Sun disk 
centers was evaluated as a function of local eclipse 
circumstances and the ratio between the lunar and 
solar radii. The estimation allows the calculation 
of the variations in the solar irradiance entering 
the terrestrial atmosphere, where the surface irra- 
diance time-patterns are simulated using a radia- 
tive transfer model. It was found that the pro- 
cedure provides a satisfactory assessment of the 
evolutionary patterns characterizing the UV solar 
irradiance during the partial eclipse of March 29, 
2006 computed by the TUV model applied to the 
Bologna (Italy) site. 
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